Kondo metal of fermionic atoms in anisotropic triangular optical lattice 
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Quantum phase transition of fermionic atoms in anisotropic triangular optical lattice was inves- 
tigated by dynamical cluster approximation combining with continuous time quantum Monte Carlo 
algorithm. The temperature-interaction phase diagram for different hoping terms and the compe- 
tition between the anisotropic parameter and interaction is presented. Our results show that the 
system undergoes Mott transition from Fermi liquid to Mott insulator while the repulsive inter- 
action reach the critical value. The Kondo metal characterized by Kondo peak in the density of 
states is found in this systems and the pseudogap are suppressed at low temperature due to the 
Kondo effect. We also propose a feasible experiment protocol to observe these phenomenon in the 
anisotropic triangular optical lattice with the cold atoms, in which the hoping terms can be varied 
by the lattice depth and the atomic interaction can be adjusted via Feshbach resonance. 

PACS numbers: 03.75.Ss, 37.10.Jk, 32.80.Hd, 67.85.Lm 
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I. INTRODUCTION 

As an ideal toolbox to investigate the quantum many- 
body phenomenon encountered in condensed matter 
physics, cold atoms in optical lattices with perfect struc- 
ture and highly controllableparameters have attracted 
great attentions recently [l|-[4j- Many cold atoms in op- 
tical lattice based quantum simulation work have been 
done with fruitful results, such as the observations of the 
Mott insulator and superfluidity of strongly interacting 
fermions in optical lattice and the quantum phase transi- 
tion between these novel stases, the quantum simulations 
of the antiferromagnetic spin chains in optical lattice, 
the frustrated Classical Magnetism in triangular optical 
lattice, Bose-Einstein condensates in optical lattices so 
forth pl-KUj. Optical lattice with different structure and 
dimension can be set up by overlapping several beams of 
laser with the same frequency and different propagating 
directions [l5l-[20j. In terms of the setting-up of the cold 
atoms in the triangular optical lattice which is geomet- 
rically equilateral but anisotropic in the hoping terms, 
three pairs of laser with the same frequency and different 
intensity in the same plane are used to create periodic po- 
tentials, as shown in Fig. 1. The fermionic 40 K atoms in 
the two hyperfine manifold are loaded into the lattice po- 
tential by combining different cooling methods including 
laser cooling and evaporative cooling. The interactions 
between fermions can be tuned by Feshbach resonance 
[Io| and the hoping between the nearest neighbor sites 
are closely related to the confining potential depth in 
different directions. 

To our knowledge, the novel states including the spin 
liquid state and the spin ice which found in the frus- 
trated systems play pivotal role in the understanding of 
the mechanism of the high temperature superconductiv- 
ity and magnetic monopoles. Hence, many experimen- 
tal and theoretical studies have focused on geometrically 



frustrated systems [21-26]. However, there is few work 
on the anisotropic triangular systems which is equilat- 
eral in the lattice constant but anisotropic in the hoping 
terms in different directions. 



In this paper, we investigate the fermionic cold atoms 
in a anisotropic triangular optical lattice by combin- 
ing the dynamical cluster approximation (DCA) [27j 
with the continuous time quantum monte carlo method. 
Even though the dynamical mean-field theory (DMFT) 
is proved to be exact in the infinite-dimensional limit and 
a good approximation even for three dimensions [28| , the 
DCA has some advantages over than the DMFT to study 
frustrated systems, because the DCA can reasonably take 
into account the nonlocal correlations of the frustrated 
systems which is overlooked in the DMFT. In the dynam- 
ical cluster approximation method, the lattice problem 
is mapped into a self-consistently embedded finite-sized 
cluster and the irreducible quantities of the embedded 
cluster are used as an approximation for the correspond- 
ing lattice quantities. This paper is organized as follows. 
In Sec. II, the artificial anisotropic triangular lattice 
with femionic atoms, the algorithm and method which 
are used to describe such artificial systems are discussed 
in detail. In Sec. Ill, we present the double occupancy 
and the density of state of the system. The Fermi sur- 
face evolution, the temperature-interaction (T-U) phase 
diagram and the competition between the anisotropic pa- 
rameters A and interaction are also shown in Sec. III. In 
Sec. V, We design a feasible experiments for the ob- 
servation of the quantum phase transitions of the cold 
fermionic atoms in the anisotropic triangular optical lat- 
tice that is geometrically equilateral but anisotropic in 
hoping terms. A brief discussion and conclusion is pre- 
sented in the last section. 
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(c)tl=t2>t3 (d)tl=t2=t3 



FIG. 1: (Color online) (a): Sketch map of experimental setup 
to form anisotropic triangular optical lattice, the degree of 
thickness of directed arrow represent the intensity of the laser 
beams and each arrow depicts the propagating laser beam 
and its direction, the circle in center of the figure represent 
the fermions. (b): cubic sketch of potentials V(x,y) in space 
for the case of isotropic triangular lattice, (c): plane sketch 
of potentials distribution of anisotropic triangular while £3<1 
and ti=t 2 =l. (d): plane sketch of potentials distribution of 
anisotropic triangular while ti=t2=ts=l. 



II. ANISOTROPIC TRIANGULAR OPTICAL 
LATTICE 

Artificial anisotropic triangular optical lattice which 
is geometrically equilateral in lattice constant but 
anisotropic in the hoping terms of different directions 
can be created through the interference of three pairs of 
laser with the same wavelength A = 1064 nm and differ- 
ent intensity in the same plane. The included angel be- 
tween each two standing wave laser should be 2tt/3. The 
fermionic 40 K atoms in the two hyperfine manifold are 
loaded into the triangular optical lattice in the same per- 
centage through evaporative cooling at a magnetic bias 
field of 203.06 G [5|. The hoping terms between nearest 
neighbor sites in different directions are closely related to 
the combining potentials depth Vo which adjusted by the 
lasers' intensity. The anisotropy in this equilateral trian- 
gular lattice can be realized through the manipulation of 
the hoping term tl, t2 and t3 in the different direction, 
as shown in Fig. 2(a). Both the non-interacting regime 
with a scattering length of a = (0 ± 10)ao and the re- 
pulsive on-site interactions between fermionic atoms can 
be tuned by Feshbach resonances, where ao is the Bohr 
radius. The confining potential in space is 

V(x,y) = Vb{[ticos(a) + 1 2 cos(/3) +£ 3 cos( 7 )] 2 

+ [tisin(a) + t 2 sin(/3) -f £ 3 sin( 7 )] 2 }, (1) 
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FIG. 2: (Color online) (a): Sketch of anisotropic triangular 
lattice, tl, t2 and t3 represents the hoping term of different 
directions respectively (b): Coarse graining procedure in the 
first Brillouin Zone when N c = 3. The red part shows the 
first Brillouin Zone. The 3 black hexagonals show the region 
of divided Brillouin Zone for DCA calculation. 



where Vo is potential depth in the x-y plane and is 
given in unit of recoil energy is E r =h 2 / (2mA 2 ) , here h 
is Planck's constant, m is the atomic mass and A is the 
wavelength of the laser beam, tl, t2 and t3 are the hop- 
ing term in different direction respectively. k x and k y are 
the two components of wave vector k = 2tt/X along x and 

y directions. a=(^§^ + ^), /3=^f^ - ^, *y=k y . 

The fermionic cold atoms with repulsive on-site inter- 
action in anisotropic triangular optical lattice can be de- 
scribed by Hubbard model and the Hamilitonian is 

H = - tijctvCj* + U Y1 n n n Hi ( 2 ) 

<ijj>a i 

where cf a and Q a represent the creation and the an- 
nihilation operator of fermios on lattice site respec- 
tively, riicr = c~l a Ci a is density operator of fermions. 
The general formula for the hoping terms is t = 
(4/v^)^r(^o/^r) 3/4 e^(-2(F /^r) 1/2 ). The different 
hoping terms in different directions can be achieved by 
adjusti ng th e lattice depth Vo in different directions. 
U = ^/S/ 7rka s E r (Vo/E r ) 3 / 4 is the on-site repulsive in- 
teraction which can be tuned by Feshbach resonance. 

In dynamical cluster approximation the lattice prob- 
lem is mapped into a self-consistently embedded finite- 
sized cluster and the reciprocal space of the lattice con- 
taining N points is divided into finite cells. The coarse- 
graining Green's function G is achieved by averaging 
Green's function G within each cell. The coarse graining 
procedure of anisotropic triangular lattice in the DCA is 
illustrated as follows. As Fig. 2 shows that the lattice 
problem of anisotropic triangular lattice is mapped onto 
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FIG. 3: The flow diagram of continuous time quantum monte 
carlo iteration loop. FT represents fourier transformation. 



three sites cluster, the Brillouin zone is divided into 3 
cells and each cell is represented by a cluster momentum 
K. The coarse-grained Green's function is 



Nc 



E 



1 



lUJ n 
k 



£ a (K,icj n ) ' 



(3) 



where summation over k is taken within the coarse- 
graining cell, the uo n is the Matsubara frequency. 

After mapping the Hubbard model onto finite size clus- 
ter, we use the continuous time quantum monte carlo 
method (CTQMC) as impurity solver to deal with the 
cluster problem. The CTQMC does not use the Trotter 
decomposition, so it is more exact than the general QMC. 
We did 10 7 sweeps in our CTQMC step. The process 
from G(K, iuj n ) = ^Ek 1 / ( iuj n ~ £ K+ £ ~ s <x(K, iu n )) 
to E c (K,io; n ) = (5(K, zcj n ) -1 — G(K,zcj n ) _1 have been 
repeated until S C (K, iuo n ) converges to required accuracy, 
as shown in Fig. 3. 



III. PHASE DIAGRAM 

We calculate the double occupancy and the density of 
state which are important to the study of the quantum 
phase transition in optical lattice systems. We also show 
the evolution of the Fermi surface's evolution at fixed 
temperature and repulsive interactions and the phase di- 
agram of the system. 

The double occupancy D occ = dF/dU = \ Xli( n *t n ^) 
as a function of interaction U for various temperature 
and different hoping terms has been investigated, where 
F is the free energy. Here we define anisotropic param- 
eter A = ts/U (i = 1,2), which describes the degree of 
anisotropy of triangular optical lattice. The double oc- 
cupancy (Docc) increases with the increase of the A = 
ts/U a t fixed interaction and temperature. With the in- 
crease of interaction the effect of temperature on Docc is 




FIG. 4: (Color online) The double occupancy D occ as a func- 
tion of the interaction U for different temperature and A has 
been shown, where A —t^/ti^ and £i=£2- (a): The double oc- 
cupancy against interaction for fixed temperature at different 
anisotropic parameter A. (b): The double occupancy against 
interaction for fixed anisotropic parameter A at different tem- 
perature. 



weakened due to the suppressing of the itinerancy of the 
atoms and suppressed while U/U is larger than 10.0 and 
12.0 for A < 1 and A > 1 respectively. As Fig. 4 shows 
that the Docc of the systems don't overlap in the fermi 
liquid phase. 

The density of state p(uj) (DOS) are calculated by max- 
imum entropy method [29j. At T=0.5, the formation 
process of a pseudogap for £3=0.8 and £3=1. 5 is differ- 
ent. For £3=0.8, a pseudogap directly formed by splitting 
Fermi- liquid- like peak at U/ti=1.0 without the appear- 
ance of the Kondo peak, while for £3=1.5, a pseudogap 
appears following after the appearance of the Kondo peak 
at ?7/£i=10.0. Still at T=0.5, the system transferred 
to insulator state while U c /U = 10-0 an d U c /U=12.0 
for £3=0.8 and £3=1.5 respectively. At T=0.25, Kondo 
resonance peak which characterized by a sharp quasi- 
particle peak with two shoulders appears at Z7 c /£i=6.5 
and U c /U=8.5 respectively at £3=0.8 and £3=1.5. A gap 
is opened when the interaction is stronger than the crit- 
ical interaction U c /ti=9.5 for £3=0.8 and U c /ti=ll.S for 
£3=1.5. The Kondo peak formed at low temperature be- 





FIG. 5: (Color online) The density of states as a function of 
frequency uj for different interaction XJ/U (i=l,2). (a): The 
density of states against frequency uj for T=0.5 and £3=0.8. 
(b): The density of states against frequency uj for T=0.5 and 
£3=1.5. (c): The density of states against frequency uj for 
T=0.25 and £3=0.8. (d): The density of states against fre- 
quency uj for T=0.25 and £3=1.5. 



cause of the effect between the local atom and the itin- 
erant atom instead of that formed at higher temperature 
through splitting the Fermi-liquid like into two parts. 

The Fermi surface's evolution for different hoping 
terms £3 at fixed £i=£2=l is shown in Fig. 6 while 
T/£ i= 0.25 and /7/£;=8.0 (i=l, 2). The spectral function 
is A(k;uj = 0) ~ — ^lim Wn ^o I mG(k,iuj n ). A linear ex- 
trapolation of the first two Matsubara frequencies is used 
to estimate the self-energy to zero frequency [30|. Due to 
the derealization of the particles the amplitude of the 
spectral weight becomes larger and the ring breadth be- 
comes narrow and distinct with the increase of the hoping 
terms for fixed temperature and interaction. The Fermi 
surface becomes a nearly fiat plane when £3 is lower than 
1.0 at T=0.25. The transition from Fermi liquid to Mott 
insulator is much faster at low temperature than that at 
relatively high temperature. 

The temperature-interaction (T-U) phase diagram and 
the competition between anisotropic parameters A and 
interaction of cold atoms in anisotropic triangular opti- 
cal lattice are shown in Fig. 7. With the increase of 
the hoping terms those corresponding transition points 
entirely shifted to larger interaction and lower tempera- 
ture. For t3=0.8, there is no Kondo peak between Fermi 
liquid and pseudogap until the T is lower than 0.45. With 
the increase of the hoping term £3, at fixed £i=£2=l, 
Kondo peak not only appears but the gap between Kondo 
peak pseudogap becomes nearly zero. When the interac- 
tion is stronger than the critical value of interaction for 
Mott transition, the system translates into Mott insu- 
lator regime which confirmed by an opened gap. The 
pseudogap are suppressed at low temperature due to the 
Kondo effect. The competition between anisotropic pa- 



FIG. 6: (Color online) The Fermi surface as a function of 
momentum k for the different hoping terms £3=0.5, 0.8, 1.0 
at fixed £i=£ 2 =l, T/£*=0.25 and U/U=S.O (i=l, 2). 



rameter A and interaction at T/£^=0.5 in Fig. 7 (d) shows 
that with the increase of the interaction,the transition 
points A for the fermi liquid to Mott insulator increase 
and the slope decreases, respectively. 



IV. EXPERIMENTAL PROTOCOL 

We propose an experiment scheme for the observation 
of the quantum phase transition of fermionic cold atoms 
in geometrically frustrated anisotropic triangular optical 
lattice. First, we use three pairs of laser in the same plane 
at wavelength A = 1064 nm with different intensity to 
create optical lattice. The included angle between each 
two beams of laser equals 2tt/3 A balanced mix- 

ture of the fermionic 40 K atoms in two hyperfine mani- 
fold states, |mF> = |-9/2> and |-7/2> states, is prepared 
th roug h evaporative cooling at a magnetic bias field 
0, [32| ^ Then a spin mixture of atoms in the |-9/2> and 
|-5/2> states can be loaded into the optical lattice poten- 
tials through tuning the scatte ring le ngth to the desired 
value by Feshbach resonance |33l-[37j. We can achieve 
the zero interaction regime with a scattering length of 
a=(0=L10)ao, where ao is Bohr radius. The repulsive on- 
site interaction U/U between fermionic atoms can be ob- 
tained by transferring the atoms in \-7/2> state to the 
|-5/2> state during the evaporation cooling. With the in- 
crease of the interactions the system will transform from 
Fermi liquid (a s <48ao) to Mott insulator(a s >lllao) at 
T=5.96nK for the lattice depth Vo=10E r . a s indicates 
the s-wave scattering length which is used to determine 
the effective interaction region. In order to detect the 
double occupancy, further tunnelling should be prevented 
through increasing the potential depth Vq of optical lat- 
tice rapidly and the energy of atoms on doubly occupied 
sites need to be shifted by approaching Feshbach reso- 
nance. One spin component of fermions on the double 
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FIG. 7: (Color online) (a): The temperature-interaction 
phase diagrams of fermionic cold atoms in anisotropic tri- 
angular optical lattice for £3=0.8. (b): The temperature- 
interaction phase diagrams of fermionic cold atoms in 
anisotropic triangular optical lattice for £3 = 1.0. (c): The 
temperature-interaction phase diagrams of fermionic cold 
atoms in anisotropic triangular optical lattice for £3 =1.5. (d): 
The competition between \=t$/ti (i=l, 2) and interaction U. 

occupied sites is transferred to an unpopulated magnetic 
sublevel by using a radio- frequency pulse. The double oc- 
cupancy can be deduced by the fraction of transformed 
atoms obtained through the absorption imagining [38| . 



After completely turning off the confining potential, the 
atoms expand anisotropically for several milliseconds due 
to the anisotropy of the potential depth in different di- 
rection. One can get the Fermi surface through the ab- 
sorption imagining [H, [39| . 



V. SUMMARY 

We investigate the quantum phase transitions 
of fermionic cold atoms in geometrically frustrated 
anisotropic triangular optical lattice with spin half 
fermions by dynamical cluster approximation method 
combining with the continuous-time quantum monte 
carlo algorithm in the frame of half-filled and one band 
Hubbard model. The system transform from Fermi liquid 
into Mott insulator while the repulsive on-site interaction 
reaches certain value for different temperature. We also 
found Kondo metal phase at middle regime temperature. 
The Pesudogap is totally suppressed at low temperature 
due to the Kondo effect. It is found that with the in- 
crease of £3 at fixed £1=^2 ? the transition lines shift to 
larger interaction and lower temperature. In the Mott 
insulator regime, the effect of the repulsive interaction 
on double occupancy becomes unapparent. 
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